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SUMMARY:  Using  one-dimensional  gas  flows  and 
quasi-two-dimensional  electrical  discharge  models,  we  calculated  the 
two-dimensional  gain  curve  contour  lines  for  equal  values  on 
perpendicular  electrode  planes.  We  studied  the  electrical  discharge 
area  gas  flow  speeds  and  electrical  discharge  parameters  as  well  as 
gas  constituents  and  similar  items  and  their  effects  on  the  rules  or 
patterns  of  changes  along  the  direction  of  the  gas  flow.  As  far  as 
each  gas  pressure  is  concerned,  they  had  flow  speed  ranges  which 
corresponded  to  optimum  gain  curves.  We  took  calculation  results  and 
experimental  results  and  carried  out  a  comparison.  >*•* — 
This  article  studies,  in  electrical  discharge  areas,  the 
electrical  discharge  transverse  flow  CW OOa  laser,  with  its  three 
mutually  perpendicular  directions  or  axes  for  light  axis,  electrical 
current  direction,  and  gas  flow  direction.  In  these  electrical 
discharge  areas,  at  the  same  time,  one  sees  the  occurrence  of 
molecular  (atomic)  ionization  and  compounding,  electron  state  and 
vibrational  state  electrical  excitation,  collision  relaxation  fading 
of  laser  stimulation  in  the  various  sets  of  dimensions,  as  well  as, 
light  absorption,  radiation,  oscillation  and  other  similar  processes. 
These  formed  a  type  of  complex,  non-equilibrium  flow  movement  problem. 
Past  theoretical  work  led  to  the  understanding  and  analysis  of  a  great 
deal.  However,  due  to  the  fact  that  the  models  were  overly 
simplified,  it  was  difficult  to  make  detailed  descriptions  of  gain 
peculiarities.  For  example,  T.A.  Cool  [l]  and  H.A.  Hassan  [2]  did 
work  in  which  it  was  assumed  that  the  electrical  excitation  occurs  in 
the  upper  reaches  of  the  light  cavity,  before  the  gas  flow  enters  the 
light  cavity  and  that  there  already  is  particle  number  inversion.  In 

the  light  cavity,  there  are  only  molecular  collision  processes  and 

-  » 

light  absorption,  excitation,  and  radiation  processes.  A.-J.  Deraaria 
[3]  and  others  only  solved,  under  stable  state  conditions,  for  the 
speeds  of  processes,  and  have  not  yet  bpen  able  to  obtain  the 
distribution  of  gain  coefficients  in  total  excitation  areas.  E. 
Armandillo  [4],  on  the  subject  of  parallel  plate  electrode  transverse 
flow  laser  theory  and  experimentation,  did  an  exhaustive  analysis. 
However,  in  order  to  overcome  the  difficulties  of  calculations,  he 
selected  for  use  a  type  of  current  flow  distribution  which  was  capable 
of  making  the  calculations  and  the  results  of  experiments  agree  with 
each  othe*-.  Moreover,  he  assumed  that  the  OOa  ionization  reached 


50$.  This  type  of  method  has  a  relatively  large  localization.  All 
the  above  work  is  set  up  on  the  basis  of  a  one-dimensional  analysis. 
But,  experimentation  still  indicates  that  the  gain  coefficient  number 
has  a  two-dimensional  distribution  [5,6].  This  article  presents  a 
type  of  simplified  one-dimensional  gas  flow  and  quasi-two-dimensional 
electrical  discharge  model,  using  numerical  value  calculation  methods 
to  quantify  or  semi-quantify  simulated  two-dimensional  gain 
distributions  along  flow  movement  directions  perpendicular  to 
electrode  planes  as  well  as  the  effects  of  flow  movement  parameters  on 
gain  distributions. 

I.  THEORETICAL  MODELS 

Electrical  discharge  models:  consider  an  electrode  structure  of 
a  form  in  which  a  metallic  plate  ABCD  is  the  anode,  EE'K'K  is  the 
cathode  (it  is  possible  that  these  could  be  tubes  or  plates  or 
needles) ,  and  the  perpendicular  distance  between  the  electrodes  AE=H 
(Pig.  1).  The  working  gases  enter  from  the  plane  AEKD.  Between  the 
two  electrodes,  one  adds  an  appropriate  voltage,  producing  florescent 
electrical  discharge,  forming  an  electrical  discharge  area  the  side 
area  of  which  is  AEPB  with  a  length  of  L.  Except  for  the  dark 
area  in  the  vicinity  of  the  cathode  EE'K'K,  and  the  dark  area  in  the 
vicinity  of  the  anode  ABCD,  the  other  large  segments  are  regular 
column  areas.  If  we  assume  that  the  gas  temperature,  the  flow  speed, 
the  pressure,  and  the  current  density,  as  well  as  other  similar  macro 
quantities  along  the  direction  of  the  y  axis  are  all  the  same,  but, 
that  in  the  two  directions  of  the  x  and  y  axes  there  are  still 
differences,  we  form  a  two-dimensional  graph  for  the  x-z  plane.  If  we 
again  assume  that,  in  the  gas  flow,  the  momentum,  mass,  and  energy 
transmission  process  along  the  z  direction  can  be  disregarded, 
approximately  taking  each  x-y  planar  layer  and  seeing  it  as  standing 
on  its  own,  it  is  only  necessary  to  give,  at  a  certain  altitude,  the 
current  density  and  the  electrical  field  strength,  and,  it  is  then 
possible,  for  that  height  or  altitude,  to  carry  out  one-dimensional 
gas  dynamics  calculations,  obtaining  one-dimensional  distributions  for 
the  x-y  plane  in  question,  showing  the  gas  flow  temperature, 
concentration,  and  small  signal  gain.  Again,  we  carry  out  calcu- 
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lations  for  different  altitudes  z,  and,  in  the  same  way,  then,  form 
our  quasi-two-diraensional  electrical  discharge  current  movement  model, 
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1.  Fig.1  Electrical  Discharge  Regular  Column  Area  2.  Light  Out  3* 
Gas  Flow 

Within  normal  column  areas,  the  current  density  j  is  distributed 
according  to  equation  (1),  and,  it  is  a  function  of  the  height  z: 


{Hb+  (o  —  b)z]L 

In  this  equation,  J  is  the  overall  current.  a=AB.  b=EF.  H  is  the 
distance  between  electrodes.  By  adjusting  a,  b,  H,  and  L,  it  is 
possible  to  flexibly  change  the  current  density  on  the  x-z  plane, 
using  the  appropriate  machine  configuration  and  placement. 

The  distribution  of  electron  density  %  is  approximately  given 
by  the  equation  below: 


In  the  equation,  vv  is  the  electron  transition  speed.  e  is  the 
electron  electric  charge. 

On  the  cathode  surface  and  in  the  vicinity  of  the  daric  space, 
electron  movement  speeds  are  very  small.  Effective  vibrational 
emmission  speeds  are  also  very  small.  Because  of  this,  the  effective 
electron  density  is  zero.  In  order  to  satisfy  this  one  condition,  we 
assume  that,  beginning  from  a  certain  place  in  the  vicinity  of  the 
cathode,  the  electron  density,  according  to  the  sine  function  numbers, 
gradates  toward  zero  in  the  direction  of  the  cathode. 

As  far  as  parallel  electrode  plates  are  concerned,  the  electric 
field  E  is  a  constant.  However,  because  of  the  fact  that,  along  the  x 
direction,  the  gas  flow  temperature  gradually  goes  up,  the  particle 
number  density  n  correspondingly  decreases.  Therefore,  this  E/n, 
along  the  x  direction,  gradually  goes  up.  As  far  as  needles  (tubes) 
go,  as  compared  to  plate  electrodes,  E,  along  the  x  direction,  goes 
down,  n  is  also  reduced  along  the  x  direction.  Experiments 
demonstrate  [5]  that  E/n,  in  regular  columnar  areas,  does  not 
experience  large  changes.  In  these  areas,  one  can  approximate  it 
with  constants. 


1.  Pig. 2  Energy  Level  Chart  2.  Electrical  Excitation  3*  Collision 
Relaxation  4-  Stimulated  Amplitude  of  Radiation 
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In  regular  columnar  areas,  the  interval  between  gas  molecules  as 
well  as  the  interval  between  molecules  and  electrons  shows  the 
existence  of  frequent  collisions.  Moreover,  carrying  out  energy 
transfer,  their  process  can  be  described  simply  by  the  use  of  Fig.  2. 

We  assume  that  rotation  and  flat  motion  take  place  at  equilibrium. 

Moreover,  we  ignore  ionization  and  the  disintegration  reaction  of 

under  the  effects  of  electrons.  According  to  traditional  methods, 

one  takes  C03  and  Na  molecules  and  simplifies  their  vibrational 

states  to  be  the  mutually  independent  simple  vibration  forms  I,  III 

and  N.  Each  vibrational  form  is  a  resonance  particle,  and,  using  the  259 

vibration  temperature  T,  in  order  to  describe  it,  the  particle  number 

uses  the  Boltzmann  distribution  expression: 


l\v, 

—  n\'"e  k‘> 


(3) 


n4f»  expresses  the  i  vibration  form  1th  vibration  energy  level 
particle  number  density.  v,  is  the  frequency.  Tt  is  the 
vibrational  temperature,  h  and  k  respectively  are  Planck's  and 
Boltzmann's  constants. 

On  the  basis  of  Fig.  2,  the  particle  number  formation  rates  for 
the  three  vibration  forms  appear  in  the  equations: 


(4) 


(5) 


(6) 


1 .  radiation  2.  electrical  excitation  3*  radiation  4*  electrical 
excitation  5*  electrical  excitation 
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Each 


quantity  can  be  seen  in  reference  [7]  at  appendix  B. 


II.  MATHEMATICAL  EXPRESSIONS 


In  regular  columnar  areas,  at  each  altitude  or  height,  the  gas 
flow  is  a  steady  state  one-dimensional  ideal  gas,  obeying  the 
conservation  equation  set  below: 


~(pu)^0 

(7) 

p«*L+ie._o 

dx  dx 

(8) 

Hi 

i 

ii 
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Hi 

(9) 

*  rfl  (*_2,  3,  K) 

(10) 

P-P^T 

(11  ) 

In  these  equations,  p,  «, p, K  ««  and  M  respectively  represent 

gas  density,  flow  speed,  pressure,  the  enthalpy  of  each  gram  of  gas, 
electrical  field  strength,  dynamic  speed  and  the  molecular  weight  of 
compounds.  Because,  here,  consideration  is  not  given  to  heat 
transfer,  and  there  is  also  no  laser  output,  it  follows  that 


dg 

dx 


-o. 


The  right  end  of  equation  (10)  is  expressed  by  equation  set  ( 4 ) — ( 6 ) . 
Take  equation  (9)  and  integrate  it  against  x.  After  rearranging,  one 
gets 


i 
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(12) 


T- 


K-HV-EU+EX 

3fhr(^co,+^N,)+  *  ^b'J 


In  this,  HV  is  the  vibration  energy, 


EX-^ 


The  subscripts  express  the  regular  columnar  entry  conitions.  is 

the  component  M  in  the  molecular  percentage  of  the  gas  compound,  to  = 

1  ,2,3  represents  C02, ,  N3>  ,  and  He.  Prom  equations  (7),  (8),  and 

(ll),  it  is  possible  to  derive: 


«“-!(*  ±  VlP-4  {R,M)T^ 

(13) 

P  =  Po + Po  «o  —  -y  Po  «0 

(£±  (R/M)T) 

(14) 

In  this*,  s~  (P0+po«o)/po«o0 

The  computational  procedure  is:  given  the  entry  parameters  2V  Po *pu 
and  the  regular  columnar  formation  a,b,H,L,  and  j ,  as  well  as  E 
distribution,  take  the  equation  set  ( 4 ) — ( 6 )  and  (12)  and  couple  them 
to  derive  n<  and  T.  Integrate  beginning  from  x=0  straight  through 
until  after  the  end  of  BCG-F.  From  equation  (3),  derive  T<.  .  From 
equations  (11),  (13),  and  (14),  derive  P,  u.  p0  .  Take  the  derived  »»ioo 
and  "ooi  and  substitute  them  into  equation  (15)*  It  is  then 
possible  to  derive  the  small  signal  gain  60  along  the  x  axis  for  a 
certain  altitude  value  z. 


In  the  equation,  A  =  10.6  microns.  ^ai  is  the  spontaneous 
radiation  life  (5-38  sec).  *•  is  the  collision  frequency.  Due  to 
the  fact  that  operational  pressure  is  relatively  high,  here, 
consideration  is  only  given  to  pressure  widening. 

Using  the  Rugge-Kutta  (phonetic  approximation)  method,  we  carried 
out  a  numerical  value  integration.  The  electrical  excitation  speed 
parameter  which  is  used,  the  electron  transition  speed,  and  the 
electron  temperature  are  all  selected  from  [8],  The  collision 
relaxation  speed  parameter  which  is  used  is  selected  from  [9-12]. 

III.  CALCULATION  RESULTS  AND  DISCUSSION 

The  primary  parameter  ranges  which  this  article  calculates  are:  T0 
=  293K.  -P0  =  20-50  (ton).  u*  =  30-310  m/sec.  The  gas  consti¬ 

tuents  are  [CO,]  ]:[N]  =  [0.03-0.05]  :[0.29-0.71  ]  .  The  rest  is  [He]. 

j  =  1 0—30  mAmp/sq.cm, 

E/n=( 1 .9-3 *0) x  10-“  volts. sq.cm. 

Pigs.  3-7  are  the  calculated  values  for  the  small  signal  gain  and 
flow  movement  parameters  along  the  direction  of  the  flow  movement  at 
locations  where  the  height  or  altitude  is  z=1 .5cm.  Prom  these  Pigs, 
it  is  possible  to  see  that  the  gas  flow  pressure  P  along  the  x 
direction  does  not  cnange  greatly.  The  temperature  T  and  tne  speed  u 
show  a  general  rise.  T  is  approximately  equal  to  T a,  .  And  TK  is 
generally  larger  than  Ts  .  If  one  continues  to  raise  the  gas  flow 
speed,  then,  the  calculated  values  for  maximum  temperature  rise  and 


maximum  speed  difference  go  down.  This  is  shown  in  Table  1 .  This 
explains  why  the  flow  speed  increase  is  capable  of  effectively 
eliminating  waste  energy.  It  is  capable  of  overcoming  the  lower 
energy  level  blockage  caused  b;7  the  temperature  increase  and  creating 
a  bottleneck  effect  for  the  gain  drop.  This  is  the  advantage  of  the 
flow  movement  laser  device  over  the  static  laser.  Because  of  this, 
the  gas  flow  pressure  P  is  capable  of  going  up  a  certain  amount. 

The  primary  difference  between  the  areas  in  Pig.  3  and  Pig.  4  is 
that  the  gas  flow  pressures  are  not  the  same.  In  Pig.  4,  the  various 
parameter  curves  along  the  x  direction  all  show  changes  which  are 
comparatively  steeper  than  those  in  Pig.  3*  This  explains  the  reason 
why  increases  in  the  pressure  have  relatively  large  effects  on  the 
laser  emission  dissipation  reaction  of  collisions  between  molecules. 
Pig.  6  shows  that,  when  the  relative  or  specific  electric  power  is  not 
greatly  changed,  the  pressure  increase  will  cause  the  small  signal 
gain  G0  along  the  x  direction  to  go  down  relatively  fast. 


1.  Pig.  3  Changes  in  ri,  u,  p,  g0  Along  the  x  Direction  E/n=2.2x  10” 

volt. sq.cm;  =0.05:0.27:0-68;  MN=5cm;  F°  =20  ton 

2.  ra/sec  3-  ton  4-  m  5>  cm 
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1.  pig.  6  Go-P0~x  Graph  E,/n=2.2x  10  volt. sq.cm.;  «* 

=  30m/ sec ;  ^  =293*;  *:*»:*‘  =1:2:7;  MN=5cm  2.  m  3-  ton 

4.  ton  5-  ton  6.  ton  7.  cm 
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1  .  Pig.  7 

=50  ton; 


Go-u<,~z 
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1 .  Table  1  The  Relationships  Between  Flow  Speed,  Rises  in 
Temperature,  and  Increases  in  Speed  2.  ton  3-  m/sec  4.  m/sec  5-  cm 
6.  cm  7.  volt. sq.cm. 


The  difference  between  Fig.  4  and  Fig.  5  lies  in  the  fact  that 
the  changes  in  E/n  along  the  x  direction  are  not  the  same.  In  Fig.  4, 
E/n  is  a  constant.  In  Fig.  5,  E/n  along  the  x  axis  goes  up  (E  is 
approximately  equal  to  a  constant.  n  along  the  x  axis  goes  down. 

This  is  with  the  use  of  simulated  parallel  plate  electrodes.)  In  Fig. 
4,  the  peak  gain  value  is  1  .0  m  .  In  Fig.  5>  the  peak  gain  value  is 
1.3  in'-  .  Also,  u  and  T  go  up  quickly. 

When  the  externally  added  electric  power  does  not  change,  from 
Fig.  7,  we  can  see  that  rises  from  30  m/sec  to  70  m/sec.  The  peak 

gain  value  is  basically  unchanged.  It  is  just  that  the  location  of 
the  peak  value  is  shifted  back.  The  effective  gain  surface  area  is 
expanded.  Because  of  this,  we  arrive  at  the  conclusion  that,  for  a 
certain  operating  pressure,  there  exists  an  optimum  flow  speed  range. 
If  one  excessively  increases  the  flow  speed,  there  is  no  attendant 
advantage.  In  Fig.  7,  P<>  =  50  tons.  The  G0~uo~x 

relationship  for  «<,  =  60-70  m/sec  is  basically  in  line  with  the 

experimental  results  in  [6].  The  difference  lies  in  the  fact  that  the 
experimental  value  for  q0  along  x  goes  down  relatively  fast.  This 
is  due  to  the  fact  that,  in  our  calculations ,  we  selected  for  use  a 
uniform  distribution  along  the  x  direction  for  current  density. 


G®o~* 


Fig.  8  gives,  for  the  same  operating  conditions,  the  G<>~x 
graph  when  there  are  different  E/n-  G0  increases  with  E/n,  and 
rapidly  rises.  However,  after  it  reaches  a  peak  value,  it  goes  down 
at  even  a  somewhat  steeper  rate.  The  amount  of  mass  flow  does  not 
change.  The  rise  in  the  value  of  E/n  signifies  an  increase  in 
specific  or  comparative  power.  At  the  current  experimental  levels,  an 
increase  in  the  specific  electric  power  can  easily  produce  electrical 
arcing,  destroying  normal  operation.  If  it  is  possible  to  overcome 
this  difficulty,  an  appropriate  increase  in  the  value  of  E/n  is 
advantageous  for  an  increase  in  gain. 


1  .  Fig.. 8 
50  tons; 


G<,~E/»~x  Relationship  Graph  >v  =  70  m/sec; 

To-  =  293K;  =  1:2:7;  MN=  >cm;  E/n  = 


(  2 .0-3 .0)  x  10- 


volt. sq.cm.  2.  m  3*  cm 
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( l)  (a) 

AB- 5.5**;  SF-2.9M*;  B- 3** 


1 .  (a)  Equal  Gain  Distribution 
Areas  AB=5*5cm;  EF=2.9cm;  H=3cm  2. 


Graph  in  Gradient 
cm  3 •  cm 


Flourescence 


1 .  (b)  Equal  Gain  Distribution  Graph  in  Gradient  Form 
Flouresence  Areas  A£=8cm;  EF=3cm;  H=3cm  2 .  era  3 •  cm 
Areas 


U 


°-2  ,0.3 


.0.5 


.78 


0.78, 


(*) 


1.  Fig.  10  Comparison  of  Calculations  and  Experiments  “»■  =70  m/sec;  -P, 

=20  tons;  r0  =29  3K;  E/n=2.2x  ,10'ie  volts,  s^.  cm, ; 

=5:27:68,  MN=5*5cm  2.  m  3*  Calculated  Value  4*  Experimental  Value 
5 .  cm 


Fig.  9  shows  the  equal  gain  curve  graphs  for  the  x-z  plane.  (a), 

(b),  and  (c)  represent  the  calculated  results  for  three  types  of  262 

different  regular  columnar  area  forms.  These  graphs  are  similar  to 
the  experimental  curves  in  [5,6].  However,  the  curve  dispersal 
intervals  are  different.  The  reason  for  this  is  that,  in  reference 
[6],  the  lower  part  of  the  cathode  has  a  section  of  cold  gas  flow 
through  it,  mixing  with  the  laser  activation  medium  and  reducing  the 
electron  density.  Because  of  this,  the  lower  area  experimental  curves 
are  drawn  together  and  are  not  terribly  similar  to  this  model.  In 
reference  [5],  because  the  pressure  is  very  high,  the  curves  are 
naturally  compressed  together. 

Fig.  10  gives  a  comparison  of  the  y=14  curves  calculated  in  this 
article  and  Fig.  4  in  [6].  The  theoretical  and  empirical  curves  for 
peak  gain  values  do  not  differ  much.  However,  the  positions  for  the 
former  peak  values  are  moved  considerably  back.  Moreover,  the  gain 
changes  behind  the  peak  value  are  relatively  even  and  slow.  The  main 
problem  is  that  this  article  assumes,  in  regular  columnar  areas,  that 
the  current  density  along  the  x  direction  is  uniform.  In  reality,  due 
to  electron  diffusion  and  incoming  flow  shock,  the  fluorescence  area 
is  moved  back,  and  along  the  x  direction  is  certainly  not  uniform. 


IV.  CONCLUSION 


Due  to  the  fact  that  we  selected  for  use  a  uniform  current 
distribution  model,  all  the  physical  quantities  created  show 
relatively  even  values  in  their  changes  along  the  x  direction.  If  it 
is  possible  to  select  for  use  even  more  appropriate  current 
distributions  and  field  strength  distributions,  this  will  cause  the 
calculation  results  to  be  even  better. 

From  calculation  results,  one  can  see  that  it  is  only  necessary 
that,  for  the  various  z  layers,  the  entry  parameters  be  the  same. 
Because  the  flow  movement  parameters  given  rise  to  by  the  introduction 
of  current  show  very  small  changes  along  the  z  direction,  taking  the 
calculation  results  at  the  position  x=4cm  as  an  example,  when  z 
changes  from  0  to  H=3cm,  and,  3  m/sec.  4T~8K  K . ,  and  AP~o  ■ 

This  explains  why  ignoring  the  momentum,  mass,  and  energy  transfer 
for  different  z  planes  is  feasible. 
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